
replace g_lat = g_lat*60*60
replace g_lon = g_lon*60*60

gen g_lat_help = (g_lat - 169028.66)/10000
gen g_lon_help = (g_lon - 26782.5)/10000

#delimit ;
gen EM =  2600072.37
         + 211455.93 * g_lon_help
		 -  10938.51 * g_lon_help * g_lat_help
		 -      0.36 * g_lon_help * g_lat_help^2
		 -     44.54 * g_lon_help^3  ;
#delimit cr

gen GKODX = EM - 2000000

#delimit ;
gen NM =  1200147.07
         + 308807.95 * g_lat_help
		 +   3745.25 * g_lon_help^2
		 +     76.63 * g_lat_help^2
		 -    194.56 * g_lon_help^2 * g_lat_help
		 +    119.79 * g_lat_help^3  ;
#delimit cr

gen GKODY = NM - 1000000
